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ABSTRACT 

We present observations of O VI A1032 emission around the helium white 
dwarf KPD 0005+5106 obtained with the Far Ultraviolet Spectroscopic Explorer. 
Previously published data, reprocessed with an updated version of the calibration 
pipeline, are included along with new observations. The recent upward revision 
of the white dwarf's effective temperature to 200,000 K has motivated us to re- 
analyze all the data. We compare observations with photoionization models and 
find that the density of the O VI nebula is about 10 cm~'^, and that the stellar 
flux must be attenuated by about 90% by the time it impinges on the inner face 
of the nebula. We infer that this attenuation is due to circumstellar material 
ejected by KPD 0005+5106 earlier in its evolution. 

Subject headings: circumstellar matter — planetary nebulae: general — stars: 
individual (KPD 0005+5106) — ultraviolet: ISM — white dwarfs 
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Introduction 



The DO (helium-rich) white dwarf KPD 0005+5106 is the h ottest known A vhite dwarf 



star. Its effective temperature was measur ed to be 120,0 K 



recently revised upward to 200,000 K by IWerner et al.l (20071) b ased on the detection of 



Ne VIII lines in its spectrum. Subsequently, IWerner et al. 



by 



Werner et al.l (Il994l ). and 



(l2008a[ ) identified emission from 



Ca X in its UV spectrum, the highest ionization stage of any element observed in a stellar 
photosphere. Weak X-ray emission at 1 keV detected with the ROSAT PSPC has been inter- 
preted as eviden ce for a hot ("^ 2 x 10^ K), fast, perhaps shock-heated wind, t hough it may 



be pho t ospheric (lO'Dwver et al.ll2003l : IChu et al. 



2004: Drake fc Werner 2005; Werner et al. 



2008bh . IChu et all J2004h find that the star is surrounded by an optical nebula, visible in 



[O III] A5007 emission, roughly 3° in diameter. They find that the nebula contains about 
70 M0 of material, and argue that this large mass implies that the material must be in- 
terstellar rather than circumstellar. Closer in, lOtte et al.l (|2004j ) discovered bright, diffuse 
O VI A1032 emission extending at least 3.5' from the star. 



The lOtte et al.l (120041 ) study is based on data from 15 lines of sight observed with the 
Far Ultraviolet Spectroscopic Explorer (FUSE) between 2001 and 2003. KPD 0005+5106 was 
observed as a calibration target throughout the FUSE mission, and spectra of half a dozen 
additional sight lines are now available. Meanwhile, the FUSE data-reduction software has 
evolved considerably, particularly with the release of CalFUSE v3. The revised effective tem- 
perature of the white dwarf, the new data, and the improved calibration pipeline motivated 
us to re-examine the nature of the O VI emitting nebula around KPD 0005+5106. 



2. Observations, Data Reduction, and Analysis 



FUSE consists of four coaligned spectrographs. Two employ optics coated with Al+LiF 
and record spectra over the wavelength range 990-1187 A; the other two use SiC coatings 
and are sensitive to wavelengths below the Lyman limit. The four channels overlap between 
990 and 1070 A. Each channel has three apertures that simultaneously sample different parts 
of the sky. The low resolution (LWRS) aperture is 30" x 30" in size. The medium resolution 
(MDRS) aperture is 4" x 20" and lies about 3.5' from the LWRS aperture. The 1.25" x 20" 
high-resolution (HIRS) aperture lies midway between the MDRS and LWRS apertures. A 
fourth location, the reference point (RFPT), is offset from the HIRS aperture by about 1'. 
When a star is placed at the reference p oint, all three apert ures sample the backgroun d sky. 
For a complete description of FUSE, see IMoos et al.l (I2OOOI ) and lSahnow et al.l (I2OOOI ). 



To identify observations for this study, we searched the Multimission Archive at Space 
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Telescope (MAST) for FUSE observations of KPD 0005+5106 obtained through the MDRS, 
HIRS, or RFPT apertures, for which the LWRS aperture should sample the background sky. 
Because its sensitivity at 1032 A is more than twice that of any other channel, we consider 
only data from the LiF lA channel. We use data obtained in time-tag mode, which preserves 
arrival time and pulse-height information for each photon event. Our final sample, consisting 
of data obtained along 21 lines of sight, is presented in Table [U the location of each sight 
line is plotted in Figure [H Our sight lines now range over more than 180° around the star. 

The data were reduced using an implen aentation of CalFU SE v3.2, the final version of 
the FUSE data-reduction soft ware package (IDixon et al.l 120071 ). optimized for faint, diffuse 
emission. Details are given in iDixon et al. J2006h . but we highlight a few points here: the 
program operates on one exposure at a time, producing a flux- and wavelength-calibrated 
spectrum from each. Background subtraction is turned off. Zero-point errors in the wave- 
length scale are corrected by fitting a synthetic emission feature to the Ly/5 airglow line. 
The spectra from individual exposures are shifted to a heliocentric wavelength scale and 
combined. Finally, the spectra are binned by 14 pixels (53 km s~^), a value chosen so that 
two binned pixels just span a diffuse emission feature. 

To improve the statistical significance of our results, we combine spectra with the lowest 
signal-to-noise ratios with those of adjacent sight lines, as noted in Table [2l The resulting 
spectra are fit with a synthetic O VI A1032 emission feature using the non-linear curve-fitting 
program SPECFIT (jKrisslll994l ). The model consists of a top-hat shaped function 106 km 
s~^ in width, representing the LWRS aperture, convolved with a Gaussian, representing the 
instrument line-spread function (~ 25 km s~^) convolved with the intrinsic emission-line 
profile. Free parameters in the fit are the level and slope of the background (assumed linear) 
and the amplitude, central wavelength, and FWHM of the Gaussian. 

Our model fits yield line intensities with I/a{I) > 2.9 for 14 spectra. The best-fit 
parameters are presented in Table [2l (Note that the sight lines are ordered anticlockwise, 
with those at radius 1.75' followed by those at radius 3.5'.) Measured intensities range from 
7.3 to 22.1 kLU. (1 LU = 1 photon s~^ cm^^ sr^^.) A fit to the sum of all 14 spectra yields an 
intensity of 15.5±0.8 kLU, a line width of 30±10 km s~^, and a local standard of rest velocity 
fLSR = 11 ± 2 km s~^. Individual lines of sight show considerable scatter about the mean 
velocity. Several sight lines have FWHM values < 25 km s~^, indicating that their O VI 
features are unresolved. Two spectra show no evidence of O VI A1032 emission; for them, 
we compute 3a upper limits to the O VI intensity by first modeling the linear continuum, 
then increasing the strength of a model O VI A1 032 line un til rises by 9 (corresponding 
to a 3cr deviation for one interesting parameter; lAvnilll976l ). Plots of selected spectra and 
best-fit models are presented in Figure [2l The asymmetric shapes of some model emission 



-4- 



features are artifacts of our sparse sampling of the model lines. 



There are two caveats. First, these observations are relatively short. To maximize the 
signal-to-noise ratio of our spectra, we use all of the data, i ncluding that obtai ned during 
orbital day. To test for contamination by geocoronal emission (jShelton et al.ll2007l ). we repeat 
the analysis, using only data obtained during orbital night. We find that, in all cases for 
which the night-only spectrum has significant exposure time (> 1500 s), the day-plus-night 
and night-only line intensities agree within the night-only uncertainties. Second, we do not 
attempt to fit the fainter component of the O VI doublet at 1038 A. It is only half as bright 
as the 1032 A hue ( in an optically-thin medium) and may be blended with interstellar C II* 
A1037 or geocoronal O I A1039. 

The reddening towar d KPD 0005-1-510 6, E{B — V) = 0.06 (§3.1), reduces the observed 
O VI intensity by ~55% (iFitzpatrickl Il999l ) . The line intensities presented in Figure [Hand 
Table [2] are not corrected f or this extinction. T he stellar spectrum shows strong absorption 
from molecular hydrogen (IMurphy et al.ll2004l ). The H2 absorption feature most likely to 
overlap the O VI A1032 emission line is the Lyman (6,0) -R(4) transition at 1032.35 A. Its 
velocity, wlsr = — 13 km s"^ (J. Kruk, private communication), places it some 100 km s~^ 
from the velocity of the summed O VI emission features derived above. Assuming that both 
the O VI emitting and H2 absorbing gas are photoexcited and have low turbulent velocities, 
then their intrinsic line widths should be about 10 km s~^ (a value consistent with the 30 km 
s~^ width of the summed O VI feature, which includes the instrumental line spread function). 
The O VI emission and H2 absorption features are thus well separated. 



Our results generally agree with those of lOtte et al.l (120041 ). Where they differ, our 
intensities tend to be lower. Exa mples are sight lines S4044403 and M1070224 (map keys 
S43 and M24), which lOtte et al.l fit separately, obtaining intensities of 25 ± 4 and 25 ± 7 
kLU, respectively. Because the intensity of the M24 feature is poorly constrained in our fits, 
we combine the two spectra and obtain a best-fit intensity of 17.0 ± 2.5 kLU. For sight line 
M1070223, the di fference is m ore extreme: our intensity of 9.4 ±3.0 kLU is half of the 18 ±8 
kLU reported bv lOtte et al.l . though their value is poorly constrained. Note that, for most 



lines of sight, lOtte et al.l use only data from the nighttime portion of the orbit. 



3. Modeling the Nebular Emission 

A star with an effective temperature of 200,000 K is expected to produce a significant 
amount of high energy radiation, capable of ionizing oxygen up to O^"*". Thus, the O VI 
emission we detect around KPD 0005+5106 is most likely from gas photoionized by the 
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star. We explore the conditions in the nebula by running photoionization models. The 
calculations, described in d e tail in §3.2, were performed with version 08.00 of Cloudy, last 
described by iFerland et al. I (Il998l ) . 



We require the models to reproduce two observational result s: the "typical" O VI in- 
tensity and the spatial extent of the [O III] emission detected by IChu et al.l ( 2004 ). As we 
will show, these conditions place strong constraints on the incident flux and on the density 
of the emitting gas. We do not seek to explain the distribution of the O VI emission, nor 
the variations in its intensity (Figured]), using these model calculations. First, we discuss 
the stellar parameters of KPD 0005+5106, including its distance. The latter is necessary for 
obtaining the spatial scales of the emitting regions. 



3.1. Stellar Parameters 



In their initial analysis of KPD 0005+5106, IWerner et al.l (119941 ) derive an effective 
temperature Tes = 120,000 K and a surface gravity logg = 7.0. Comparing these parameters 
with post-asymptotic giant branch (post-AGB) evolutionary tracks, the y estimate the star's 



mass to be M* = 0.59 Mq. Comparing its V-band flux (V = 13.32; iDownes et al.l Il985l ) 
with model predictions, they derive a distance to the star of 270 pc. At this distance, the 
3.5' radius of the outer set of FUSE sight lines (Figured]) corresponds to a physical extent 
of 0.27 pc. 

The detection of Ne VIII absorption features in the F USE spectrurn of KP D 0005+5106 
requires an upward revision of its effective temperature. I Werner et al.l (120071 ) find that the 
Ne VIII lines are best fit by models with T^s = 180,000 - 200,000 K and log^ = 6.5 - 7.0. 
Subsequently, Ca X emission lines were identified in the F USE spectrum, and fi ts to these 
lines yield T^s = 200,000 - 220,000 K and log^ = 6.2 - 6.5 JWerner et a l.l l2008ah. For both 



specie s, the best-fit value of the surface gravity rises with effective temperature. IWerner et al. 



( 120071 ) point out that the star's He II lines, in both the UV and the optical, are better fit 
by a model with T^s = 200,000 K and log^ = 6.5 than by the original T^s = 120,000 K, 
logs' = 7.0 model. 

K. Werner has kindly provided us with a synthetic spectrum of the star. It is based on a 
Teflf = 200,000 K and logg = 6.5 model atm osphere composed of He, C, O and Ne, with C = 
0.003, O = 0.0006, and Ne = 0.01 by mass (IWerner et al.ll2008ah . J. Kruk has provided the 
fully-reduced FUSE spectrum of KPD 0005+5106. From these, we can estimate both the 
star's interstellar reddening and it s distance. We be gin by scaling the synthetic spectrum 
according to the reddening law of iFitzpatrickl (119991 ) for various values of the color excess 
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E(B — V). We scale the result to match the star's V-band flux. Finally, we compare the 
predicted and observed flux at 1060 A in the FUSE spectrum. We flnd that the FUSE 
data are best fit by E{B ~ V) = 0.06, a value similar to the E{B — V) = 0.07 derived by 
Downes et al.l (119851 ) from lUE data. 



As IWerner et al.l (l2008bl ) point out, an effective temperature T^s = 200,000 K and a 
surface gravity log q = 6.5 place KPD 0005+5106 on the 0.7 Mq post-AGB evolutionary 



track of 
radius R 



Wood fc Faulknerl (119861 ). Adopting this mass and surface gravity, we derive a 



0.08 Rq. The scale factor required to make the reddened model (at 5500 A) 
match the star's V-band fiux is SF = 1.4 x 10~^^. Give n the star's radius and this scale 
factor, its distance is simply = tiR^/SF (lKurucg|l979l ). or d = 820 pc. At this distance. 



the 1.75' radius of the inner set of FUSE sight lines corresponds to ~ 0.42 pc (1.3 x 10 cm), 
the 3.5' radius of the outer set corresponds to ~ 0.84 pc (2.6 x 10^^ cm), and the 1.5° radius 
of the [O III] nebula corresponds to ~ 21 pc (6.6 x 10^^ cm). 

We repeat the calculations using a second synthetic spectrum (also provided by K. 
Werner) with T^fi = 180,000 K and logg = 6.5. Scaling the model to match the observed 
V-band and 1060A fiuxes, we derive a color excess E{B — V) = 0.06, which is the same as 
for the hotter model. Assuming the same stellar mass as before, we derive a distance to the 
star of 780 pc. 

Small changes in the surface gravity affect the stellar radius, but do not significantly alter 
the overall shape of the emergent continuum. A star with Tgs = 200,000 K and logg = 6.2 
would have a radius of 0.11 -R0, which is about 40% larger than that for logg = 6.5. The 
luminosity would therefore be greater by about 90% (since L oc R^), and the derived distance 
to the star would increase to 1.2 kpc. 



3.2. Cloudy Models 

Cloudy is a versatile code that can be used to calculate the ionization structure and 
emission from photoionized gas under a broad range of conditions. Our models for the O VI 
nebula detected by FUSE are relatively simple. We assume an ionization-bounded slab of 
uniform density illuminated by the white dwarf. The main input parameters are the shape 
of the ionizing continuum, the incident flux on the cloud surface, the density structure of the 
gas, and the elemental abundances in the gas. Calculations are performed in one dimension 
assuming spherical symmetry; however, the effective geometry can be made plane parallel 
by choosing a large radius of curvature. 

In our models, the abundances of elements are held flxed. We use the default solar 
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values as specifi ed in Cloudy: C and O from AUende Prieto et al. (2001, 2 002); N, Ne, Mg, 
Si and Fe from iHolwegerl (1200 ll ): and the rest from iGrevesse fc Sauval I (Il998l ). We take 
the Tcff = 200,000 K synthetic spectrum of KPD 0005+5106 as the shape of the incident 
continuum. (In §3.2.2 we use the Teg = 180,000 K spectrum to study how the results depend 
on the effective temperature of the star.) The incident flux and the gas density are allowed 
to vary (§3.2.3). 

The physical distance from the star to the nebula scales with the distance from the star 
to the earth. The stellar distance is derived by scaling the model spectrum to match the 
observed fluxes in the V-band and at 1060 A. Therefore the flux normalization at the nebula 
does not depend directly on the distance between us and the star, but rather on the stellar 
flux observed at earth and the angular separation between the star and the inner face of the 
nebula, which we take to be 1.7'. 



3.2.1. Initial model 

We flrst explore the case of the unattenuated stellar continuum incident on the O VI 
emitting nebula. The inner face of the nebula is assumed to be 1.7' from the star, placing 
it just within the inner set of FUSE observations (clearly it can be no further away). This 
corresponds to a distance of 1.3 x 10^^ cm at 820 pc. The normalization of the continuum may 
be described in terms of the flux of photons with energy greater than or equal to 113.9 eV, 
which is the energy required to ionize O^^ to 0^+. (This ionization energy corresponds to a 
photon frequency of about 2.76 x 10^^ Hz.) From the synthetic spectrum of KPD 0005+5106, 
we flnd that the luminosity of such high-energy photons is ~ 2.8 X 10^6 photons s'^ . The 
flux of these photons at the inner face of the nebula is ~ 1.4 x 10^ photons s~^ cm~^. For 
this initial model, we choose a hydrogen number density of 1 cm"'^, which is a typical value 
for the interstellar medium (ISM). 

The calculations predict that the O VI emitting zone extends about 20 kpc, and that the 
O VI A1032 emissivity is about 1.5 x 10~^^ erg s~^ cm~^ over most of the region. O VI A1032 
is a resonance line, and the Cloudy predictions of emissivity take into account radiative 
transfer effects via the escape probability formalism. Although the escape probabilities are 
calculated for a one-dimensional slab, we assume that the predicted emissivity holds for 
our observed lines of sight, which lie perpendicular to the direction of the stellar flux vector 
through the model nebula. Furthermore, we assume that the path length through the nebula 
is about 4.2 x 10^^ cm, which is roughly equal to the extent of the O VI emission in the plane 
of the sky (about 5.7', the distance between sight lines M08 and M23 in Figured]). Using 
these values, we flnd that the expected intensity from the model nebula is ~ 5 x 10~^ erg 
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cm~^ sr~^. A color excess E{B — V) = 0.06 (§ 3.1) toward KPD 0005+5106 corresponds 
to an extinction at 1032A of about a factor of 2 (iFitzpatrickl Il999l ) . Therefore, the model 
predicts that the observed intensity would be about 2.5 x 10~^ erg s~^ cm~^ sr~^, which 
equals about 13,000 LU. This value is consistent with the FUSE observations. 

Although this model reproduces the observed O VI intensity, it does not correctly re- 
produce the extent of the optical nebula. In the model, the [O III] emission arises beyond 
the O VI, at a distance of between 20 and 32 kpc from the star. This is about t hree orders 
of ma gnitude larger than the ~ 20 pc radius of the [O III] nebula detected by IChu et al. 



3.2.2. Teff = 180,000 K model 

The lower limit on the effective temperature of KPD 0005+5106 is 180,000 K. We 
calculate a model using the Teg = 180,000 K synthetic spectrum to examine the sensitivity 
of the nebular properties to the stellar effective temperature. All other model parameters 
remain unchanged from the initial model. 

The luminosity of 0^+ ionizing photons is ~ 1.4 x 10^^ photons s~^, about half that of 
the 200,000 K stellar spectrum. Their flux is ~ 6.8 x 10^ photons s~^ cm~^ at the inner face 
of the nebula. 

The O VI zone extends about 12 kpc, just over half the distance predicted by the 
200,000 K model. The total O VI flux from the 180,000 K model nebula is only half that 
of the 200,000 K model. However, the intensity in the inner few kpc is about 20% higher. 
This is because more of the 0^+ ions are ionized to O^"*" by the 200,000 K spectrum than 
the 180,000 K spectrum. The [O III] emission arises at a distance of between 12 and 22 kpc 
from the star. This is still almost three orders of magnitude larger than the observed radius 
of the optical nebula. 

If the effective temperature of KPD 0005+5106 were at the lower limit allowed by the 
Ne VIII lines, the unattenuated flux impinging on gas at typical interstellar densities would 
still produce an optical nebula much larger than observed. If the temperature were at the 
higher limit (220,000 K), the discrepancy would be increased. In the rest of this paper, we 
consider only the best-fit 200,000 K synthetic spectrum for KPD 0005+5106. 
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3.2.3. Fiducial model 

We ran a grid of models using the T^s = 200,000 K synthetic spectrum for the continuum 
shape, varying the incident flux level and the density of the emitting gas. The incident flux 
levels ranged from the unattenuated stellar continuum used in the initial model down to 
one-tenth this value. (The latter corresponds to 1.4 x 10^ 0^+ ionizing photons s~^ cm^^.) 
The hydrogen number densities ranged from 1 cm~'^ to 100 cm~^. 

For models in which the continuum is unattenuated, a gas density riu of about 40 cm^'^ 
is required to produce an [O III] nebula with a radius of about 20 pc. In that case, however, 
the predicted O VI intensity is over 200 times the observed value. As the incident flux is 
lowered, the density required to reproduce the observed radius of the [O III] nebula also 
decreases. For each such combination, the predicted O VI intensity decreases with incident 
flux and density, matching the observed values only when the incident flux is ~ 10% of the 
unattenuated value. The corresponding hydrogen density in the gas is 10 cm~^. We take 
this to be the fiducial model for the nebula around KPD 0005+5106. 

In Figure [3] we plot the predicted emissivities of O VI, [O III], and Ha against depth into 
the cloud. The average O VI emissivity out to a few parsec is about 2 x 10~^^ erg s~^ cm~'^. 
This emissivity is slightly higher than that predicted by the initial model and, under the 
same assumptions used above, yields an observed intensity of 17,300 LU. The O VI nebula 
would be visible out to a few parsecs, i.e. beyond the outer ring of FUSE observations. The 
[O III] emissivity falls rapidly beyond about 24 pc, and the Ha emissivity extends slightly 
beyond 26 pc . The total [O I II] to Ha flux ratio is about 6. This is consistent with the 



conclusion of IChu et al.l (120041 ) that the Ha nebulosity associated with the nebula around 



KPD 0005+5106 is faint compared with the [O III]. 

The presence of photoionized O VI emission around a hot white dwarf is not unexpected. 
As shown by our model calculations, a 200,000 K star can maintain oxygen in this highly 
ionized state out to a distance of several kiloparsecs for typical ISM gas d ensities and can pro 



duce the observed O VI A1032 intensity. (A similar point was made by lPanagia fc Terzian 



19841 . who showed that UV radiation from hot white dwarfs can maintain the presence of 
C^"*" in the Galactic halo.) The puzzling result is that, in order to produce the observed 
O VI emission and an [O III] emitting zone confined to a radius of ~ 20 pc, the ionizing 
continuum must be attenuated by about 90% between its emergence from the stellar surface 
and its incidence on the inner face of the nebula. In the next section, we discuss how this 
attenuation may be explained by the presence of circumstellar material around the white 
dwarf. 
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4. Discussion 

4.1. Circumstellar material around KPD 0005+5106 

The prior evolution of KPD 0005+5106 - and in particular its mass-loss history - are 
unknown. However, the star is expected to be surrounded by material lost while it was on 
the asymptotic giant branch. We make the assumption that a shell of material surrounds the 
star and examine the effect of such a shell on the stellar continuum. To do so, we calculate 
Cloudy models with the KPD 0005+5106 {T^s = 200,000 K) spectrum impinging on a shell 
and compare the emergent continuum with the incident. 



From the lWood fc Faulkneii (119861 ) evolutionary track, we estimate that KPD 0005+5106 
left the asymptotic giant branch some 2000 years ago, after ejecting its outer envelope. We 
thus fix the inner radius of the shell to be 1 x 10^^ cm, the distance traveled in 2000 years 
by material moving at a typical PN expansion velocity of 20 km s~^. Assuming that enve- 
lope ejection took a few thousand years, we fix the outer radius of the shell at three times 
the inner radius. The density, assumed to be uniform, is treated as a variable parameter. 
Abundances are fixed at the solar values. We find that the continuum above 113.9 eV is 
attenuated by about 90% when the density of the gas in the shell is uh = 3500 cm~'^. The 
total mass in the shell (assuming that it covers An steradians) is then ~ 0.4 Mq, comparable 
to the 0.2 Mq typical of planetary nebulae. 

The shape of the emergent spectrum differs from that of the KPD 0005+5106 synthetic 
spectrum in several ways. There is a sharp decrease at the He II edge at 54.4 eV from which 
the spectrum begins to recover at ~ 100 eV. While the O^"^ ionizing continuum is reduced 
by 90%, the continuum between 54.4 eV and 100 eV is reduced by over 99%. The continuum 
between 13.6 eV and 54.4 eV is reduced by only about 70%. 

We next calculate a Cloudy model of the nebula using this emergent spectrum propa- 
gated out to the inner face of the nebula (1.3 x 10^^ cm from the star). As described above, 
the shell density is chosen such that the flux of O^"^ ionizing photons incident on the nebula 
is the same as that for the fiducial model. The hydrogen density in the nebula is 10 cm~'^, 
also equal to that in the fiducial model. The difference in the shape of the continuum has 
an impact on the predicted properties of the photoionized nebula. 

The resulting model has an O VI zone that extends out to only about 2 pc, which is 
much smaller than that predicted by the fiducial model (Figure [3]), but still extends beyond 
the outer ring of Ft/S'i? observations (Figured]). The O VI intensity in the inner 1 pc (which 
encompasses the FUSE sight lines), however, is about 3 times that of the fiducial model. 
The [O III] zone extends out to 45 pc, about twice that predicted by the fiducial model (and 
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thus about twice the observed extent of the optical nebula). 

The O VI intensity and the extent of the [O III] zone can be brought closer to the 
observed values if the flux emergent from the shell is reduced by an additional 32%, such 
that there are ~ 9.5 x 10'' O^^ ionizing photons s^^ cm~^ impinging on the nebula. A model 
with this lower flux yields an O VI zone extending about 2 pc from the star (as above). The 
O VI intensity in the inner 1 pc is within 5% of that predicted by the fiducial model, and 
the [O III] zone extends out to about 28 pc. 

The gas density in the shell, its inner and outer radii, and the abundances of elements 
are all unconstrained parameters of the models. Additionally, the nebular gas density may 
be adjusted. Because there are no independent observational constraints on the putative 
circumstellar shell around KPD 0005+5106, more detailed calculations are not warranted. 

We have shown (in §3.2) that in order to explain the O VI and [O III] observations of 
the nebula around KPD 0005+5106, the stellar flux must be attenuated. We have shown 
that this attenuation can be achieved by a shell of gas with properties consistent with mass 
loss from the star while on the asymptotic giant branch. We have further shown that the 
emergent spectrum is capable of photoionizing a nebula and producing the observed O VI 
and [O III] emission. However, we do not make any attempt to model both the shell and the 
nebular emission self-consistently. A shell such as the one we have assumed would produce 
its own emission. To our knowledge, there has not been a careful examination of the region 
around the star at sufficiently high resolution to place strong limits on any such emission. 

4.2. Spatial distribution of the O VI emission 

It is clear from Figure [1] that the O VI emission has an azimuthal structure, with the 
faintest emission from the southern and northwestern edges of our sample and the brightest 
towards the southwest. The pattern is as clearly seen in Tabled where the intensities exhibit 
a monotonic increase followed by a monotonic decrease going around the inner and outer 
rings. The only exception is sight line M34. The M34 value may be less reliable than the 
others because it has the shortest exposure time and is barely a 3a detection. If we ignore 
M34, then all the sight lines outside of the cone between lines connecting the white dwarf and 
M09 on one side and M24 on the other have O VI intensities less than or equal to 9.5 kLU. 
Within this cone, all sight lines have O VI intensities greater than or equal to 14.5 kLU. At 
the edges of the cone, only the inner ring shows the brighter O VI; M08 and M23 are fainter 
(Figure [1]). 

The observed distribution of the O VI intensities could be due to an azimuthal depen- 
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dence of either (i) the gas density in the nebula or (ii) the incident ionizing flux. We examine 
the conditions necessary for each of these alternatives, restricting ourselves to the predictions 
from one-dimensional models with uniform gas density. 

For a given incident flux, the ionization parameter (ratio of ionizing photon density to 
gas density) decreases with increasing density, while the emission measure increases with 
increasing density. Thus, there is some value of density at which the O VI emissivity peaks. 
Our flducial model assumes a gas density of 10 cm~^. Raising the density to 15 cm~^ 
reduces the O VI emissivity by a factor of 1.5. If the density is 20 cm~^, then the emissivity 
falls by a factor of 2.2. These emissivities yield observed intensities (under the path length 
assumptions used in §3.2) of 11.5 kLU and 7.9 kLU, respectively. Higher-density models also 
predict smaller and brighter [O III] nebulae than the flducial model doe s, but no obvious 



azimuthal structure is seen in the [O III] images presented by IChu et al.l (120041 ). A model 
with a density of 1 cm~'^ predicts an O VI emissivity that is half the value of the flducial 
model, but is ruled out because its [O III] zone extends to 2 kpc. 

If the gas density is flxed, then a lower incident flux will yield fainter O VI emission. 
An incident flux 90% that of the flducial model results in a factor of 1.4 drop in the O VI 
intensity. If the incident flux is 80% that of the flducial model, then the intensity drops 
by a factor of just over 2. These values correspond to observed intensities of 12.5 kLU and 
8.5 kLU respectively. The [O III] intensity and the extent of the [O III] emitting region are 
not strongly affected by changes in the incident flux - they change by factors of less than 
20% among these models. 

Spatial variations in the flux incident upon the O VI emitting nebula could be achieved 
by moderate density variations in the circumstellar shell around KPD 0005+5106. For 
example, in the shell considered in §4.1, if the density were 4000 cm~^ rather than 3500 cm~^, 
then the 0^+ ionizing continuum would be attenuated by about 98%, rather than 90%. The 
circumstellar shell is not an opaque object casting a shadow, but a source of diffuse radiation, 
with considerable angular extent when viewed from any point on the inner face of the nebula. 
Therefore we cannot simply infer the azimuthal structure of this shell based on the observed 
distribution of O VI intensities. Calculating the radiative transfer through the circumstellar 
shell to derive the incident continuum at each point in the nebula is beyond the scope of this 
work. 

Of the two alternatives for explaining the observed azimuthal distribution of O VI emis- 
sion, the latter, i.e., ionizing flux variations, is the one that has less effect on the distribution 
of [O III] emission. Therefore, we suggest that it is the more straightforward explanation. 

One implicit assumption in all of our models is that the O VI and [O III] arise in a 
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single nebula with uniform density. This need not be the case. The mass loss of stars on 
the asymptotic giant branch is a complex process; in the case of KPD 0005+5106, up to 
2 Mq may have been lost over a period of 10^ — 10^ years prior to the final mass ejection 
that resulted in the circumstellar shell we have inferred in §4.1. If such were the case, then 
it would be necessary to consider (at least) an outer region of circumstellar material and the 
interstellar medium as independent emitting regions. One may then naively expect that all 
the O VI arises in the former, while the [O III] emission is dominated by the latter. Without 
further data, detailed explorations of such a scenario are premature. 



5. Concluding Remarks 

In this paper, we have presented a new analysis of the O VI A1032 emission around the 



hot h elium-rich white dwarf KPD 0005+5106. At the time of our earlier analysis (jOtte et al. 



20041 ). the star was believed to have an effective temperature of 120,000 K and li e at a distance 



of 27 pc. Since then, the effective temperature of the star has been revised (IWerner et al. 



20071 ). and consequently its inferred distance has increased to 820 pc (§3.1). 



We have rao deled our O VI measurements and the [O III] A5007 emission around the 



star ( jChu et al.l 120041 ) as a photoionized nebula. We have used a synthetic spectrum of 
KPD 0005+5106 for the ionizing continuum and treated the nebula itself as a uniform, 
one-dimensional slab. Based on the comparison of models with data, we have inferred 
the existence of a circumstellar shell between the star and the O VI emitting nebula that 
attenuates the stellar flux by about 90%. We find that the nebula has a density of ~ 10 cm~'^. 

The distribution of O VI intensities on the plane of the sky shows a clear azimuthal 
structure. This structure may be due to a variation either of the gas density or of the 
incident flux. Within the framework of our simple models, we present the conditions under 
which we may expect the observed O VI intensity distribution for each alternative. We 
find that variations in the density would affect the brightness and the extent of the [O III] 
nebula in ways that probably would be observable, but which are not seen. However, a more 
careful examination of the optical emission with accurate photometry is required before any 
definitive statement can be made. The O VI intensity is more sensitive to variations in the 
incident flux, and the observed differences can be explained without affecting the properties 
of the predicted [O III] emission. 

Our results clearly highlight the need for a thorough examination of the region around 
KPD 0005+5106. One possible tracer of high ionization gas is [Ne V] A3426. Our Cloudy 
models predict that the strength of this line is about 10 times that of O VI A1032. A 
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photometric, narrowband imaging study of the nebula in [O III] and [Ne V] emission would 
yield diagnostic information about the distribution of material and the ionization structure 
of the nebula. We have inferred the existence of a circumstellar shell based on the synthetic 
spectrum of KPD 0005+5106 and photoionization models for the nebular emission. Its 
existence is plausible considering the evolutionary history of stars on the asymptotic giant 
branch. A direct observation of such a shell would be of great value. 
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X Offset (arcmin) 

Fig. 1. — Observations of diffuse O VI emission around the hot white dwarf KPD 0005+5106. 
Colors represent the measured O VI A1032 intensity in units of kLU (1 LU = 1 pho- 
ton s^^ cm~^ sr~^) recorded through the FUSE LWRS aperture. For M17 and M26, the 
colors represent 3a upper limits. The radius of the outer ring is 3.5 arcmin. Labels refer to 
the map keys defined in Table [H North is up and east is to the left. 
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Fig. 2. — Selected spectra showing the region around O VI A1032. The data are binned by 
14 pixels, and best-fit model and error spectra are overplotted. For M1070217, the model 
spectrum represents a 3cr upper limit. Interstellar C II * A1037 and O VI A1038 are present 
in some spectra, as are the geocoronal [O I] AA1028, 1039 lines. 
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Fig. 3. — Emissivities of selected lines predicted by the fiducial model for the nebula. The 

O VI emissivity translates to an observed intensity of 17,300 LU. consistent with the FUSE 
data. The [O III] emission extends out to about 27 pc, consistent with optical observations. 
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Tablc 1. Observations 



Observation 






Position Angle** 


Distance 


Map 


ID 


Date 


(degrees) 


(arcmin) 


Key 


iVllU 1 UZiUU 


2001 


Sep 29 




3 ^ 


A/ro'^ 


A/T1 n7n9nfi 

iVllU 1 u^uu 


2001 


Sep 29 


9fi9 


1 R 
1 . o 


A/TOfi 

iVlUU 


M^ 070908 


2002 


Sep 23 


ouu 


o. o 


A/IOS 

iVlUO 


Ml 070900 


2002 


Sep 23 


9QQ 


X . o 


A/TOQ 


A/Tl 07091 1 


2002 


Oct 21 




o.o 


Ml 1 


Ml 07091 9 


2002 


Oct 22 




1 S 

l.O 


Ml 9 


Ml 07091 A 


2002 


Oct 23 








Ml 07091 


2002 


Oct 23 




1 8 


Ml ^ 


M1070217 


2002 


Dec 17 


168 


3.5 


M17 


M1070221 


2003 


Sep 07 


311 


1.8 


M21 


M1070223 


2003 


Dec 03 


181 


3.5 


M23 


M1070224 


2003 


Dec 03 


181 


1.8 


M24 


M1070226 


2004 


Jul 24 


321 


3.5 


M26 


M1070227 


2004 


Jul 24 


321 


1.8 


M27 


M1070229 


2006 


Sep 02 


272 


3.5 


M29 


M1070234 


2006 


Dec 26 


151 


1.8 


M34 


S4053801 


2001 


Oct 23 


256 


2.2 


SOI 


S4053802 


2004 


Nov 23 


215 


2.2 


S02 


S4053803 


2003 


Dec 02 


207 


2.2 


S03 


S4053804 


2004 


Jul 23 


346 


2.2 


S04 


S4054403 


2002 


Dec 22 


188 


2.2 


S43 



'^Position angle of LWRS aperture, measured east of north, 
relative to the star. 
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Table 2. Measured O VI Intensities and Upper Limits 



Sight 


Time*^ 


Intensity'' 








FWHM'^ 


Line 






Q /l\Tc 

of r\ 


(km s~ 


) 


(km s~ ) 


M34 


1836 


18.8 ±6.2 


3.0 


-13 ± 


23 


120 ± 60 


S43/M24 


6147 


17.0 ±2.5 


6.9 


4± 


7 


10 ± 30 


S03 


9596 


19.3 ±2.7 


7.3 


8± 


7 


50 ± 30 


S02 


14800 


22.1 ±2.0 


11.1 


18 ± 


5 


50 ± 20 


M12/M15 


8029 


19.4 ±2.8 


7.0 


11 ± 


7 


50 ± 40 


sm 




1 7 n -1- 1 8 




7 + 


K 

o 


7 + zLD 


M06 


3581 


16.1 ±3.8 


4.2 


13 ± 


11 


20 ± 80 


M09 


4330 


14.5 ±3.5 


4.1 


8± 


16 


120 ± 50 


M21/M27 


6048 


7.6 ±2.3 


3.3 


40 ± 


14 


10 ± 110 


S04 


All! 


8.7 ±2.2 


4.0 


4± 


11 


20 ± 110 


M17 


3160 


< 9.0 










M23 


2521 


9.4 ±3.0 


3.1 


8± 


16 


10 ± 140 


M11/M14 


6346 


17.8 ±2.5 


7.2 


8± 


7 


10 ± 30 


M05/M29 


4457 


16.7 ±2.8 


5.9 


4± 


8 


10 ± 100 


M08 


4214 


7.3 ±2.5 


2.9 


26 ± 


9 


3± 30 


M26 


2806 


< 9.5 











'^Total exposure time. 

^1 LU = 1 photon s~^ cm~^ sr~^. 

^Statistical significance of quoted intensity, I/a{I). 

'^Gaussian FWHM values include the smoothing imparted by 
the instrument optics. Values less than ~ 25 km s~^ indicate 
that the emission feature is unresolved. 



